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Abstract 

In this work we solve exactly a class of three-body propagators for 
the most general quadratic interactions in the coordinates, for arbi- 
trary masses and couplings. This is done both for the constant as the 
time-dependent couplings and masses, by using the Feynman path in- 
tegral formalism. Finally the energy spectrum and the eigenfunctions 
are recovered from the propagators. 
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1 Introduction 



In a very recent work, A. Chouchaoui [1] calculated the propagator for the 

problem of three identical particles in one dimension. This was done by 
working in the Feynman path integral formalism [2] [3] , so expanding the class 
of three-body systems solved through this formalism as, for instance, those 
discussed by Khandekhar and collaborators [4], Govaerts [5], and others. 

On the other hand, the interest in solving problems involving time-dependent 
systems has attracted the attention of physicists since a long time. This 
happens due to its applicability for the understanding of many problems in 
quantum optics, quantum chemistry and others areas of physics [6]-[17]. In 
particular we can cite the case of the electromagnetic field intensities in a 
Fabry-Perot cavity [6] . In fact this kind of problem still represents a line of 
investigation which attract the interest of physicists [16]-[20]. 

Here we intend to expand the class of exactly solvable path integral prob- 
lems, by including the case of a general three-body quadratic interaction in 
the coodinates, both in the case of constant as in the case of time-dependent 
couphngs. This is going to be done through a suitable combination of the 
Jacobi coordinates with a further decoupling one. 



2 A general exactly solvable quadratic inter- 
actions for three-body with constant cou- 
plings 

The model which we are going to treat is represented by the following La- 
grangian 

+ cTi (fa - fi) ■ (fa - fi) + 0-2 (r2 - fi) ■ (fg - fs) + (1) 
+(73 (fa - fi) • (fa - f2) + gi ■ (f2 - fi) + §2 ■ - n) + 93 ■ {rs - ^2) ] • 

As far as we know, this problem was considered only for particular cases, 
in general without the crossed terms proportional to the a coupling constants 
[21], and most part of time for one-dimensional identical particles under 
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isotropic harmonic interactions [1],[4], [5]. As a matter of fact, we should say 
that, in fact the parameters gi, g2 and can really be time-dependent as 
we are going to see below in the text. 

For the first step, in order to reach our goal, we perform a change of 
coordinates system going to the so called Jacobi coordinates, which are char- 
acterized in this case by the following set of transformations [21]: 



( xA 



b mi 
mi2 



bm2 
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where a, b are arbitrary constants, used in order to put the Lagrangian into 
a more convenient form, and M = nii + m2 + rris is the total mass of the 
system. Note that for the sake of path-integral evaluation, we must compute 
the Jacobian of the transformation due to the impact of the transformation 
over the measure of the path- integral, and in for this transformation it is given 
simply by {ab)^. As observed in [21], this transformation is not sufficient to 
decouple the system, which will looks like 
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with the resulting masses and coupling constants given by 
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with mi2 = mi + m2- 

Note that the center of mass coordinate has decoupled entirely from the 
other ones, so that one can now worry only about the remaining coordinate 
variables. Now one can take two possible routes. The first, which is the one 
usually taken, sometimes with errors as observed in [21], by particularizing 
the problem through a constraint over the parameters: A = 0. Another 
possible and more general route, is the one we are going to choose here. In 
such route we perform a further transformation, as proposed in [22] and used 
also in [21], keeping in mind the need of decoupling the final system through 
the elimination of the crossed term. This is reached by a simple additional 
dilation and rotation transformation [22], which in the present situation is 
defined as 




where m is an arbitrary parameter with dimensions of mass, and (f) is the 
rotation angle which are going to be fixed in order to disentangle the system. 

After using these transformations, one achieves the following Lagrangian 
for the last two variables 



dt 



J^m I dYj 



dt 



aY^ + (3Yi + -fYi-Y2- 



+ F^{t)-Yi + F2{t)-Y2\, 
with the transformed constant couplings being defined as 



a 



P 



m uJl cos (j)^ + muj2 sin 0^ 



Am 



sin (20); 



mul sin0^ + muji cos0^ H — sin (20) ; 



VM1M2 



7 = m {u\ — (J^ sin (20) + 



2Am 



cos (20) ; 



(6) 



(7) 



and 



Fi = 



(8) 



F. ^ Asin0+J— /,cos0. 



Furthermore, due to the above transformation we get a change of the path- 
integral measure given by a cubic power of the Jacobian of the transformation 
(J3= (^SSy). Now we are in conditions to ehminate the crossed term in 
these last coordinates, and this is done by imposing that the angle should 
obeys: 

Solving the above equation, one can easily to note that two solutions appear, 
but they just interchange the role of the new vectors Yx and Y2 in the La- 
grangian, which lead us to conclude that both conduce to the same physical 
consequences, in such a way that we only need to work with one of them. 
We will use the following solution 



COS( 



1 / ^ ^jMxM2{ujl-ulf 

2 V \lA\^ + MxM2{ujl-ujlf 



(10) 



Using the above solution, the Lagrangian will be finally set decoupled 
and given by 



m 



E 

3=^ 



dt 
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where the final decoupled frequencies can be written as 

4A 



2 I 2 



M1M2 
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(11) 



(12) 



5 



and the forces 

Flit) ^ 

where we defined that 



m 



2 Ml 
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I m 
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{1-R)f2{t) , 



(13) 
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lAX^ + MiM^iul-ulY 

Finally we are left with the task of solving the corresponding Feynman prop- 
agator for a system of uncoupled forced harmonic oscillators. For this we can 
use the well known solutions of this system [2] . 

Remembering the fact that the final Lagrangian is a direct sum of in- 
dependent ones, the corresponding propagator must be just the product of 
three independent propagators in terms of the final variables, so we get 

3 



m 

xK2{Yi%Yi;T) Ks{n\n;T), (15) 

where the prime and the double prime denotes that a function is being eval- 
uated at the initial or at the final instant of the time interval respectively. 
Besides t = tb — ta is the time interval between these instants. The first prop- 
agator K3 is that of a free three-dimensional particle with the total mass of 
the system of three particles, and can be easily obtained [2] , as 

(16) 



— * — * 

once Y3 — Xs- The other two have the same form of a three-dimensional 
driven harmonic oscillator [2]: 

3 



Kj{Yr,Yj;T) 



m Qi 
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(17) 
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where 



GAr) ^ 



dt 



Yf-F,{t)) sm[%(t-t„)] + 



/ 2 ^ 



(18) 



dt' 



ta 



dt (Fj (t) ■ Fj {t')) X 
X sin [flj ( tb - t)] sin [flj {t' - ta)] , 



and j — 1, 2. Now, by using the above expression, one can finally write the 
final expression of the quantum propagator of the three-body system under 
examination, simply multiplying the three resulting propagators, and further 
taking back the original physical variables, through 



(V^cos0 -(y5)sin 
( /mT) sin (.[M£)cos(t) 




(19) 



and, also using the initial transformation (2). Here, however, we avoid to 
write the final expression for the sake of conciseness and, because at this 
point it is only a simple task of substituting and multiplying the decoupled 
propagators. 

From the above expressions we can extract the corresponding wave func- 
tions and energies. For this we remember that the propagator can be obtained 
from the following spectral summation 



K{z",z';T,0) = J2rn{z',t')i,^{z",t"). 



(20) 



n=0 



On the other hand, we can use the Mehler's formula [25] , 



exp 



oo 



- + b')] E -iHn (a) (6) , 

n=0 



exp [- (a^ + 62 _ 2 a 6 c) / (1 - c^) 

(21) 

in order to [22] recover the corresponding wave functions. After length but 
straightforward calculations one can finally obtain a wave function which 
consists of a three-dimensional free particle with continuous energy, and three 
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three-dimensional driven harmonic oscillators, with discrete eigenenergies. 
Once more, in order to be concise, we will write below only the quantized 
part of the wave function. 
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where a = 1,2,3 = x,y, z, S = sin(0), C = cos(0), i7„.^(-) is the riia th 
Hermite polynomial, and 



Via {t) = 



[I' diF,a (0 sin {fii iC - Q) sin {Q, (4+ -(i^|) 



\^/rn sin (Qj i) 
+ J^'" dC Fia (0 Sin {Qi (e - t„)) sin (fi^ (t^ - t)) 
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Particularly in the case of non-driving forces (^i = ^2 = ^3 = 0), the eigen- 
energies are given by 



from which we can observe that, beyond the usual degeneracy which hap- 
pens when J2a=i ^la = iVi, or I]a=i ^2a = -^2, with A^i and A^2 being integer 
numbers. In the case where the final frequencies have a rational relation, fur- 
ther degeneracies will appear. As the transformed frequencies are functions 
of the original ones and of the couplings, one concludes that, when certain 
relations between the parameters hold, the system becomes more degenerate, 
signalizing the appearance of hidden symmetries. 

3 A general exactly solvable quadratic inter- 
actions for three-body with time-dependent 
couplings 

In this section, we extend our calculation in order to include the case of time- 
dependent couplings in the three-body system. In this case the number of 
works in this matter is even more scarce and, up to our knowledge, restrict 
to the one- dimensional case [4], [5]. In this new situation, as a consequence 
of the fact that if we try to do the same above transformations, additional 
terms would appear rendering the system unsolvable if the masses were time- 
dcpcndcnt and wc must keep the masses constant. Furthermore, we should 
restrict ourselves to treat the case of isotropic frequencies, and perform some 
new transformations because the decoupled system is still time- dependent 
[11]. In view of these arguments, we deal with a system characterized by the 
Lagrangian density (1) with arbitrary constant masses, perform the trans- 
formations (2) getting a Lagrangian density with the same form of (3) but 
now, one has time-dependent frequencies and couplings. However, in or- 
der to guarantee the exact solvability we must impose the following set of 
constraints 




'2 , 



(25) 



Ml = 



mim2 



= M2 = 



II, UJi {t) = UJ2 {t) , 



(26) 



aP"m-i2 
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which leads us to fix one of the arbitrary constants through the relation 



(27) 



On the other hand, the constraint among the frequencies implies into a re- 
striction over the time-dependency of one of the coupling parameters, given 

by 



LU it) = Ujf = 

mim2 



mi2 



K21 + -\- (mlK^i + mlK32) + 
mf2 ^ ^ 

aim2 - a2mi - > = (28) 

V mi2 /J J 



M 



(X31 + K32 + 2 as) 



171317112 

Leading to the following Lagrangian 



2 \ dt fr{2 \ dt 2 



(29) 



+ X{t) X^-X2+ h{t)-X, + f2{t)-X2 



Now, we can decouple the coordinates X-i and X2, by performing a 7r/4 
rotation around the third vector like 



Xi 
X2 



X2 



(30) 



and this leads us to the following Lagrangian 



^ M dX, 
2 dt 



Ejjj ( dec j \ ^(x 



dt 



VLi {tf XI + Q2 {tfxi 



2 ;S2 



(31) 



+61 {t) • f 1 + ^2 {t) ■ X2 
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where the frequencies and hnear force- type couphngs of the decoupled three- 
dimensional oscillators are respectively given by 

(tf ^ u {tf ± e, W = ^ ( /i {t) ± /; {t)) ; i = 1, 2. (32) 

Furthermore, the path integral measure is invariant under this transfor- 
mation. As a consequence, we are left to solve the problem of a free particle 
(coordinate X3), and two forced oscillators with time-dependent frequencies 
and forces [11]. Once the case of the oscillators is more involved, we review 
briefly a way to get the propagator for them. First of all, we perform a 
translation like 

yi = Xi +f]i, (33) 

and then impose the elimination of the linear terms, which make it necessary 
to restrict the arbitrary functions ffi through the equation 

and this allow us to write the Lagrangians as 

with 

^^^/^^ •(%-2%)-^/V,(0 -ffAOd^, (36) 

From the above, we conclude that the translation used has been able 
to reduce the problem to that of a harmonic oscillator with time-dependent 
frequency. At this point we can, for instance, map the problem into that of a 
free particle, what can be done by using a set of coordinate transformations, 
including a time substitution, introduced by Cheng many years ago [23], [9], 

Zi = yi sec [/x, (t)] , Ui = tan [a^ (t)] , (37) 
where ai {t) and some auxiliary variables Sj (t) should obey the equations 

Si + Q,i {tf Si = \, si Oii = 1. (38) 
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After these transformations, we ends with the following Lagrangians 



^ 2 \dui, 



d Fj {y,. ij,. t) dGj{si,t) 



dt 



dt 



(39) 



with 



G,{s^,t) = lfi^Lm[2fij{t)]-2 



Sj Uj 



(40) 



Looking to the above equations, one can conclude that we have really 
mapped the original problem into that of free particles. Now, using the 
solution of the free particle oscillator [2] and the Van Vleck-Pauli formula 
[24] , it can be shown after straightforward calculations that the ith propagator 
becomes 



^ia {vL y'ia^ 'T, 



2 7ri sin(5,„ 



exp 



2h 



Via— + Via iVia-^Via) 
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X exp 



I fx 
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X exp [-^l ^ia (0 Via (0 d^) , 



(41) 



where it was define for the sake of compactification of the expressions that, 
for any function g{t), one have that g{T) = g" and g{0) = g'. Besides, in 
the above expression 5ia = a'l^ — a'ia- 



lmim2MiM2 



iM 



M 



2 7ri T/ 

^^^^ 



Let us now recover the wave functions through the use of the decomposi- 
tion of the propagator in terms of the wave functions, as it is given from (20), 
and again using the Mehler's formula (21), with the necessary identifications 
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^-^Ui and Q = exp(—i5i), we obtain for the discrete 



part of the wave functions 

nia,n2a {-^lai -X2a) CXp 
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X exp 
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Finally we would like to remark that some other systems which are some- 
what more general than that appeared until now in the literature could 
be included in this analysis. This is the case, for instance, of a three- 
dimensional generalization of the potential considered in the interesting work 
of Chouchaoui [1] , which would be represented by the following Lagrangian 



^ 2 



1 



K21 (r2 - fi)^ X31 (fa - fi)^ X32 (fg - f2)^ 

(45) 



dt j 2 

Ox if 2 - fi) • (ra - fi) (T2 (f2 - fi) • (fs - fa) 

^1 



+0-3 (ns - ^1) • (fs - ^2) + 



<^ ^\2 ' 
J2 - Tl) _ 
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and with some restrictions over the potential parameters. Furthermore, we 
could also include in the list of the systems which can be decoupled through 
the use of the transformations above described, some cases where particular 
three-body interactions could added to the general quadratic coordinate one. 
More precisely, one could accrescent something like 

Ay = -rJ^ + .+ 



(f2 - fi) (mi2f3 - (mifi + m2r2)) {miTi + m2r2 + mafa) 

to the interacting potential. But, once again, restrictions over the parameters 
of the potential should be imposed in order to keep the exactness of the 

propagator and wave functions solutions. Moreover, the case of the general 
quadratic spatial interactions, which we considered in this work, has the 
advantage of presenting covariant interactions between the three particles 
involved. For this reason, we prefer do not treat these cases here. 
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